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Using a generalized order parameter gradient expansion within density functional theory, we 
derive a phase-field-crystal model for liquid crystals composed by apolar particles in three spatial 
dimensions. Both the translational density and the orientational direction and ordering are included 
as order parameters. Different terms involving gradients in the order parameters in the resulting 
free energy functional are compared to the macroscopic Ginzburg-Landau approach as well as to 
the hydrodynamic description for liquid crystals. Our approach provides microscopic expressions 
for all prefactors in terms of the particle interactions. Our phase-field-crystal model generalizes 
the conventional phase-field-crystal model of spherical particles to orientational degrees of freedom 
and can be used as a starting point to explore phase transitions and interfaces for various liquid- 
crystalline phases. 
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I. INTRODUCTION 

The traditional Landau theory of phase transitions in 
which the free energy is expanded in terms of a conve- 
nient order parameter predicts the order and the scaling 
behavior of bulk phase transitions [l| in mean field ap- 
proximation. The same idea can be used for spatially 
dependent order parameters in a free energy functional 
where gradient expansions lead to a Landau-type de- 
scription of equilibrium interfaces between two coexisting 
phases (see e.g. [2[)- This approach was very successful 
for liquid-gas transition and fluid-fluid phase separation 
in mixtures and was generalized to nematic and smectic 
liquid crystalline phases by de Gennes S 0] . It was also 
extended to freezing by a multiple [Sl-Q or single order 
parameter theory. The latter involves a gradient expan- 
sion up to fourth order in the density field and leads to 
the so-called phase-field-crystal (PFC) model ^^. De- 
pending on the parameter combinations, the PFC model 
leads to stable periodic density modulations and to con- 
stant densities both in two and three spatial dimensions 
|ll| . Therefore it has been used for large-scale numerical 
investigations of statics and dynamics in the crystalline 
state including: the structure and free energy of the fluid 
crystal interface [12^ 13| , crystal growth dynamics into a 
supercooled liquid [lj|, the structure I^ and dynamics 
[ig of grain boundaries, and the Asaro-Tiller-Grinfeld 
instability [T7l - l20| . The PFC model can be derived from 
microscopic density functional theory (DFT) B l2lU24J 
which describes crystallization in equilibrium |25l - [28| us- 

_a Landau expansion in terms of density modulations 
l29j . In two spatial dimensions, this derivation was 
recently generalized to liquid crystalline phase which pos- 
sess orientational order [3d| and within quite few parame- 
ters a rich topology of the equilibrium phases were found. 

In this paper, we generalize the derivation of the phase- 
flcld-crystal model for apolar orientational degrees of 



freedom to three spatial dimensions. We start from mi- 
croscopic density functional theory and perform a gra- 
dient expansion in three order parameters namely the 
translational density and the orientational direction and 
ordering (or equivalently the nematic tensor). The result 
for the static free energy functional is richer than in two 
dimensions [30[ . The prefactors of the different gradient 
terms are expressed as integrals over microscopic correla- 
tion functions allowing thus a microscopic determination 
in terms of the interparticle interactions. The results 
are compared to macroscopic approaches which provide 
a framework of possible gradient terms allowed from gen- 
eral symmetry principles. This phase-fleld-crystal model 
generalizes the phase-fleld-crystal model of spherical par- 
ticles @-[i3l to orientational degrees of freedom. It can be 
used as a starting point to explore phase transitions and 
interfaces for various liquid-crystalline phases, in partic- 
ular including plastic and full crystalline phases where 
the translational density shows a strong ordering. 

The paper is organized as follows: in Sec. II, we derive 
the PFC model from density functional theory by ex- 
panding the orientational dependence of the density field 
up to the first nontrivial order and performing a gradi- 
ent expansion in the translational coordinate. Then, in 
Sec. Ill, we discuss the relation to symmetry-based ap- 
proaches. We finally discuss possible extensions of the 
model to more complicated situations and give final con- 
clusions in Sec. IV. 



II. DERIVATION OF THE 

PHASE-FIELD-CRYSTAL MODEL FROM 

DENSITY FUNCTIONAL THEORY 

Our derivation of the PFC model uses the micro- 
scopic static density functional theory for liquid crystals. 
We consider N particles with ccnter-of-mass positions 



{Ri \i = 1, . . . , N} and two oricntational degrees of free- 
dom each. The actual orientations of these particles are 
described by a set of unit vectors {ui \i ~ 1, . . . , N}. In 
three dimensions the system has a total volume V and is 
kept at a finite temperature T in the domain V C R-^. 

A pair interaction potential W{Ri — R2tUi,U2) be- 
tween two particles 1 and 2 is assumed. We henceforth 
consider apolar and uniaxial particles. The apolarity im- 
plies the following symmetries: 



W{r,Ui,U2) = W{-r,Ui,U2) = W{r,-Ui,U2) 
= W{r,ui,-U2) = W{r,U2,Ui) ■ 



(1) 



Examples for W{Ri — i?2, ^i, ^2) comprise excluded vol- 
ume interactions (e.g. hard spherocylinders 13 ll |32| or 
hard ellipsoids [Sj] ) , Yukawa segment models [3J-[3a] and 
Gay-Berne potentials [37l - l39l |. 

The inhomogencous one-particle density 



N 



p{R,u)^ (^5(i?-i?0^("-w0 



(2) 



provides the joint probability density to find particles 
at center-of-mass-position R with orientation ii. For an 
observable A the symbol 



{A)=^j^d^Rr 



X ^exp 
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W{Ri - Rj,Ui,Uj) 
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denotes the normalized canonical average, with the clas- 
sical canonical partition function Z, the unit sphere S2 — 
{x gM.^ : \x\ = 1} and Boltzmann's constant fee- Apolar 
particles involve the symmetry p{R, u) ~ p{R, — m). 

From classical density functional theory of inhomoge- 
ncous fluids we know about the existence of an excess free 
energy functional J^3xc such that the Landau free energy 
functional 



n{T, V, p, [p{R, ii)]) = J-id(T, V, [p{R, u)]) 



+ T,UT,V,[piR,u)]) 



d^R / d'^upp{R,u) 
V Js2 



(4) 



is minimal for the equilibrium density field for a given 
chemical potential p, temperature T and domain V. The 
ideal rotator gas functional J^id is known exactly: 

J-id(r,V,[p(i?,u)]) =kBT [ d^R [ d^up{R,u) 

Jv JS2 (5) 

X (\n{k^p{R,u)) -l\ . 

Here, A is the thermal de Broglie wavelength. The ex- 
cess free energy functional Fa^c{T,V^ [p{R,u)]) incorpo- 
rates all correlations and is not known in general, but 



there are several approximations available. Well known 
approximations include the Onsager functional 



J-cxc(T,V,[p(i?,7i)]) 



^ / d^R, / d^R2 



X / d 111 d U2p{Ri,ui)piR2,u2) 

' 02 '' S2 



(6) 



X exp 



WiRi-R2,Ul,U2) 

knT 



-1 



that becomes asymptotically exact in the low density 
limit 140[, the mean-field approximation 



T,^,{T,V,[p{R,u)]) 



1 



d^Ri I d^R2 I d^iii 



(7) 



X / d U2W{Ri- R2,Ul,U2)piRl,Ul)piR2,U2) 

IS2 



which becomes asymptotically exact at high densities for 
bounded potentials |4l| and the Ramakrishnan-Yussouff 
approximation [25| 



T^UT.VMRM) 



knT 



d-'Ri / d''i?2 

V Jv 



d^U2c''^HRl-R2,Ui,U2) (8) 



X / d^ui 
' S2 "' S2 

X (p{Ri,ui)-'p](p{R2,U2)-p 



This approximation is used in the following. It can be 
viewed as a truncated density expansion in the density 
difference p{R, u) —J) around a mean density p with the 
direct correlation function c*^^' of a reference fluid. More 
accurate forms for J^oxc for hard particles are given by 
weighted-density-approximations |42. |43| or follow from 
the fundamental- measure theory |28| . As a further ap- 
proximation we only consider weak anisotropics in the 
orientations. Thus, the leading terms in the density 
parametrization for uniaxial particles are 

p{R, u) « p(i + MR) + MR) PiiMR) ■ ")) (9) 

with P2(a^) = ^{3x^ — 1) denoting the second Legendre 
polynomial. In this expression, the real-valued dimen- 
sionless orientationally averaged density is represented 

by 



^i(i?) = T^ / d^u{p{R,u)~p) 
47rp Js2 



(10) 



It is identical to the original treatment of the PFC model 
@, Q . The dimensionless field 



MR)^T^ d^upiR,u)P2{MR)-u) , (11) 
47rp Js2 

on the other hand, measures the local degree of oricnta- 
tional order. For apolar particles, the leading anisotropic 



contribution is the third term in the big brackets on the 
right-hand-side of Eq. ©. Furthermore, the unit vector 
field Uq(R) defines the local director of the orientation 
field H- 

To begin with the derivation of the static free energy 
functional J- = J-id + J-cxc we insert the truncated ex- 
pansion ([9]) into Eq. ([5]) and expand the logarithm up to 
third order. After performing the angular integration the 
approximation 



and 



Jv 



X 2^^! + i't 
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(fR (13) 



is obtained. Next, we derive the excess free energy 
functional. For this purpose we insert Eq. ^ into the 
Ramakrishnan-Yussouff approximation ([8]) and decom- 
pose the direct correlation function into spherical har- 
monics Y;_m(w) up to second order: 



.(2) 



{R,U1,U2)^ E E ^hLmun..iR) 



lj—0 nij— — lj 



(14) 



X Y/i_™i(wi)Y/2,,„2(w2) • 
The angular integration leads to the final expression 



J^CKc[i>l,1p2,Uo] 



-kBTp^Sn^ I d''i?i / d''i?2 

'V JV 



21 i 



E E 

^j— raj— — 2lj 
l^j^2 l^j^2 



^ ^2h,2l2,mi,m2^^^ 



R2 



(15) 



X Y2;i,,„i {uoiRl))llJh+liRl) 
X Y2l2.r,i2 {uoiR2)) 1pl2 + liR2) 

Here, the expansion coefficients 



JS2 JS2 

X Yi^^rm{ui)Yi^^„i,_{u2) 



(16) 



arc not independent. If they are further decomposed into 
a series of spherical harmonics with respect to the remain- 
ing orientational unit vector (i?i — R2)/\Ri — -R2I, it is 
possible to use the identities 

00 ij I 

c^'^'>{Ru,ui,U2) = ^ oji,j^j{R) ^ ^ 

^1,^2,^—0 'mj — — lj m— — l 

KK2 (17) 

X C{li,l2j,rni,m2,m)Yi^^„ii{ui) 

X y 12,7,12(^2) Y 1, 7n{u) 



Wil,i2j(-R) 



47r 



min{/i,i2} 
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J f d^iii f d^u2 Yl 

•^■5'2 "'■52 m=-min{lid2} 



X C{li,l2,l, m, -TO, 0) Yhjn{ui)Yl2~r,i{u2) 

X c(2)(i?e2,ui,'U2) 



(18) 



to derive necessary relations between the expansion co- 
efficients from the properties of the Clebsch-Gordan co- 
efficients C{li,l2,l,Tnijm2,m) [45| . Now a gradient ex- 
pansion is performed [29| up to fourth order in the ■01 V"! 
terms of Eq. (|15p and up to second order in the f/'i V'2 and 
il'2'<p2 terms. We assume that the highest gradient terms 
ensure stability. By partial integration and omission of 
surface terms one obtains the result 



1 



J^cxc[tpi,tp2,uo]~^ I d^R 



X Ui^2^^2(V^l) +v43(A^l) 
+ ^lV'i + ^(^l+2i^2)(VV^2)' 

+ i{Ki~K2){uo-yi'2y 

+ B2 (VtAi • VV'2 - 3(uo • Vi/'i) (mo • VV'2) 

- 3?/'2 V-01 • ( (uo • V) uo + uo ( V • uo) ) 
+ 2V'2VV'2 • ((i^2 - i'l) ((uo- V)uo) 

+ -(2i^i+A'2)(wo(V-wo)) 

+ ^|(^l(V-Mo)'+i^2(wo-(VxUo))' 



+ Ki{uox{yxuo)y 



for V = M"^ . The coefficients are given by 



Ai = -kBTpHnM^^l , 
A2 = fcBVfM(^), 



B. 



An 



15V5 



M, 



(2) 

0,2 I 



A-i = ksTf ^{9M^^), - 16M^'l + lOM^'^ 



(19) 



B, = -kBTf ^ (ufl - 2m(°] + 2Mg) , ^20) 

-fcBV 



and depend on the moments 
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(21) 



of the expansion coefficients of the direct correlation func- 
tion, where e-z is the cartesian unit vector codirectional 
with the ^;-axis. These moments in turn depend on the 
thermodynamic conditions expressed by T and p. 

The functional ([T9l) contains a few simpler models as 
special cases which follow by setting ipi, -02 and uq suc- 
cessively to zero or to a constant. In detail, ii ^^i =0 
and ■02 and uq are constant, the PFC model corresponds 
to a completely isotropic system and all terms in Eq. 
([T^ vanish or are constant {^ Bi). Alternatively, when 
■01 ~ 0, "02 is constant and Mn is space-dependent, the 
Frank free energy functional [46| is recovered. Up to the 
common prefactor 0| the Frank constants for splay, twist 
and bend are Ki, K2 and K3 = Ki, respectively. The 
more general case with 0i = and space-dependent 02 
and Mo can be referred to as constant-density approxima- 
tion. If then Uq is a constant unit vector, the constant- 
density excess free energy functional reduces to a simple 
gradient expansion of second order for "02, where only 
the terms four to six arc not vanishing. Otherwise, for 
a space-dependent uq, the PFC model corresponds to 
the Landau-de Gennes free energy [4g for uniaxial nc- 
matics. More complex models are recovered for non- 
constant densities, i.e. for -01 = ipi{R). If only tJji is 
space-dependent, the PFC model has no orientational 
dependence and is equivalent to the three-dimensional 
extension of the PFC-model of K. R. Elder and cowork- 
ers [23 for isotropic particles. With a space-dependent 
Mo, this model is extended to the free energy for uniax- 
ial nematics with a constant amount of ordering. Also 
for non-constant scalar fields 0i and 02 but a constant 
ncmatic director Uq, the model is much simpler than the 
full excess free energy functional in Eq. ([T^ . because the 
computationally expensive terms that describe the cou- 
plings of Vipi and VV'2 with derivatives of the nematic 
director as well as the very involved Frank free energy 
drop out. 

In the full excess free energy functional ([T^ all these 
special cases are properly comprised. This new func- 
tional clarifies the relation between already existing sim- 
pler PFC models, contains the appropriate couplings of 
the fields -01, -02 and uq, relates the constant prefactors 
of the terms in Eq. (|19p to the direct correlation function 
and is therefore the main result of this paper. 



III. RELATION TO SYMMETRY-BASED 
APPROACHES 

In this section we make contact between the three- 
dimensional PFC model for liquid crystals based on den- 
sity functional theory with two macroscopic symmetry- 
based approaches, namely the Ginzburg-Landau descrip- 
tion and generalized hydrodynamics. The goal is to com- 
pare the central results of this paper summarized in Eqs. 
([T9l) to ((2T|) with corresponding results from Ginzburg- 
Landau analysis appropriate as a mean field description 
in the vicinity of phase transitions and the hydrodynamic 
description applicable for long wavelengths (continuum 
approximation) and low frequencies. 

For the contributions associated with density varia- 
tions and their gradients given in the second line of Eq. 
([19)) this can be done very easily. From the Ginzburg- 
Landau description for the smectic A - nematic descrip- 
tion we have for the corresponding terms [47j in the en- 
ergy density 



^«IV'P + i&i|v,vP 



^62|A0p 



(22) 



where we have used the smectic order parameter -0 = 
'0oexp(— 10) with magnitude -00 and phase introduced 
by de Gennes 0, |4a|, which is directly proportional to 
density variations (compare, for example [48|). Compar- 
ing the second line of Eq. (|T^ and Eq. (|22p , we can iden- 
tify Ai, A2 and A3 with a, 61 and 62, respectively. 

In addition, we can also make contact with the bulk 
description of smectic A, where one uses the layer dis- 
placement u parallel to the layer normal as a hydrody- 
namic variable [4^,|49|. For the gradient terms associated 
with the layer displacement, u, which is proportional to 
phase changes, one has in the energy density 



25(Vf)' 



^KiVlur 



(23) 



where in Eq. ([23)) the contribution ^ i? is associated with 
the compressibility of the layering and the contribution 
^ K is connected with bending of the layering. Since the 
macroscopic description is dealing with the bulk of the 
smectic A phase, the uniaxial anisotropy is reflected in 
the terms given in Eq. ([^5)1 . 

For the terms associated exclusively with orientational 
order in Eq. ([TO)) in lines 3, 4 and 7-10 we start the com- 
parison with the continuum description of the bulk phase 
for which we have for the analogous terms in the energy 
density 

A1(V • nf + K2{n ■ [V x n]f + Kz{n x [V x n\f 

+ L||(n,V,5)2 +L^5jj(V,5)(V,5) (24) 

+ M{V.S)[5%n,+5%n,]{V,nu). 

In Eq. ([M)) the first line is connected to gradients of the 
director field, n. It contains splay, twist and bend and 
goes back to Frank's pioneering paper [4a, [SO] ■ Lines 2 



and 3 are associated with gradients of the nematic modu- 
lus, 5', and with a couphng term, ^ M, between gradients 
of the director and gradients of the modulus [5l|, [54I . 

The Frank free energy can be easily compared with 
lines 9 and 10 of Eq. ((T9)) . We identify ip2 in the last 
section with the nematic modulus, S, and Ug with the 
director field, n. The splay constant Ki in Eq. 1^^ reads 
2i^i?/;|, the twist elastic constant in Eq. ((24)) reads 2i^2'02 
and the bend elastic constant in Eq. ([M]) is 2Kiip2- A 
similar comparison can be performed for the expressions 
given in lines 3, 4, 7 and 8 of Eq. (|T^ . We obtain for the 
contribution -- M in Eq. ([24]) ■i/;2(A'2 - Ki) and for the 
gradient terms of the nematic modulus L|| and L± are 
given by L|| = |(4i^i - K2) and L^ = |(i^i + 2i^2). 

We thus arrive at the conclusion that instead of 6 in- 
dependent coefficients in the symmetry-based continuum 
description only two independent ones are left over in the 
PFC model. This reduction in the number of indepen- 
dent coefficients is a well known feature of approximate 
approaches. It is even known in fields such as superfluid 
■^He where one finds a reduction from 6 to 3 coefficients 
for analogous terms [SlJ, l53| . In this connection it turns 
out to be quite instructive to compare the terms associ- 
ated with gradients of the orientational order in Eq. (fTO)) 
with the gradient terms in the Ginzburg-Landau descrip- 
tion of the nematic - isotropic phase transition given by 
de Gennes Q with the same terms also occuring in the 
Ginzburg-Landau description of the smectic A - isotropic 
phase transition. 

In the vicinity of the isotropic - uniaxial nematic phase 
transition one takes traditionally the second order trace- 
less tensor Qij [3 as the order parameter. It vanishes 
in equilibrium, < Qij >= in the isotropic phase and 
assumes in the uniaxial nematic phase the structure [46[ 
Qij = S{niiij — 2^ij)' with d = 2 in two dimensions and 
d = 3 in three dimensions. The gradients in Qij give rise 
to the gradient energy 



Fq = 



d RLijkiran^iQjkl^lQmn) 



(25) 



order by pointing out that the first term on the third line 
in E n. ([T^l) is the analog of the term ~ /3i in Eq. (2a) of 
Ref. jlf"^ 

Next we compare the results given in Eq. ([T^ for 
the coupling terms between gradients of the density and 
gradients of the orientational order with the results of 
the two macroscopic symmetry-based approaches. These 
terms are listed in lines 5 and 6 of Eq. ([T^ and are pro- 
portional to B2. For spatial gradients in the director 
field coupling to spatial variations in the density we find 

MM 



F„^ 



d^RXPiV^pMln, + d%nA (V,n,) (27) 



where the transverse Kronecker delta projects onto the 
plane perpendicular to the preferred direction n: (5,-J = 
^ij ~ rLirij. By comparison with Eq. (jl9|) we find A'' = 
— 3i?2^2- Finally we have for the terms coupling gradi- 
ents of the order parameter modulus to gradients of the 
density [5a | 



Fsc 



d3i?iVf,(V,5)(V,p) 



(28) 



where the second rank tensor iV is of the standard uni- 



axial form iVf; 



N^TiiTij + iVjJ,-''. A comparison with 



Eq. dUl) yields iVi = -2^2 and N2 = B2. The cou- 
pling terms listed in Eqs. (j27p and ([28)) exist in both, 
two and three spatial dimensions. Thus in comparison 
to the hydrodynamic description of the bulk behavior, 
which is characterized by three independent coefficients, 
we find one independent coefficient in the PFC model. In 
the framework of a Ginzburg-Landau approach using the 
orientational order parameter Qij we find in the isotropic 
phase 



Qrho 



d3i?P«(V,0,fc)(V,p)(%4i+<5.fe<5,/) (29) 



In the uniaxial nematic phase the tensor Lijkimn has six 
independent coefficients. Using the decomposition, for 
example, in three dimensions Qij = S{ninj — -^Sij) one 
obtains in total six coefficients as above: 3 Frank coef- 
ficients for the pure deformations of the director field, 2 
coefficients for the deformations of the modulus S, and 
one coupling term between gradients of the director and 
gradients of the modulus. In the isotropic phase this ex- 
pression reduces to [J| 



Fc 



d3i?[Li(V,Q,fc)(V,Q 



Jk) 



(26) 



i2(V,Q,fe)(VjQjfc) 



and thus to the same number of independent coefficients 
as in the PFC approach given above. We close the dis- 
cussion of the terms associated purely with orientational 



and thus one independent coefficient - as has also been 
the case for the PFC model. 

In conclusion we find that the PFC model expand- 
ing the generalized energy up to quadratic order in ori- 
entational and density variations and their gradients 
can be compared easily with corresponding terms in 
the symmetry-based Ginzburg-Landau description in the 
isotropic phase. It also emerges that the number of in- 
dependent coefficients obtained in the framework of a 
PFC model based on a density functional approach typ- 
ically contains a smaller number of independent coeffi- 
cients than predicted from generalized hydrodynamics. 
This feature is shared by other approximate approaches 
such as BCS for superfluid ^He, but has the advantage 
to predict concrete values for these coefficients, which are 
left as unknown parameters in a hydrodynamic descrip- 
tion. 



IV. CONCLUSIONS AND POSSIBLE 
EXTENSIONS 

In conclusion, we derived a phase-field-crystal model 
for liquid crystals in three dimensions from density 
functional theory. Two approximations are involved: 
first the density functional is approximated by a trun- 
cated functional Taylor expansion similar in spirit to the 
Ramakrishnan-Yussouff theory. Then a generalized gra- 
dient expansion in the order parameters is performed 
which leads to a local free energy functional. There are 
three order parameters, namely the translational density 
which corresponds to the scalar phase-field variable ■01 in 
the traditional phase- field-crystal model, an orientational 
direction given by a three-dimensional unit vector uq and 
the nematic order parameter ■02- In the two latter quan- 
tities the gradient expansion is performed up to second 
order while it is done to fourth order in ipi . This ensures 
that the traditional phase-field-crystal model [^, Q is re- 
covered as a special case in which there is no orientational 
dependence of the full density. The additional terms are 
all in accordance with macro scop ic approaches based on 
symmetry considerations [52|, [53 . The full static free en- 
ergy functional allows for a wealth of stable liquid crys- 
talline phases. How the phase diagram depends on the 
prefactors should be explored by further numerical stud- 
ies. Once the stable phases are known, the structure 
of interfaces between two coexisting phases can be ad- 
dressed, not only the isotropic-nematic interface |56| - [58| 
but also interfaces which involve one or two translational 
ordered crystalline phases. 

The analysis presented here can be generalized or ex- 
tended to quite a number of different situations. First 
of all one could in principle include higher order gradi- 
ents both in ipi as well as in uq and 02- The former is 
in particularly mandatory if a more realistic description 
of the translational crystalline density field is targeted 
which is highly peaked in a three-dimensional solid [59[ 
or if higher-order orientational distributions should be 
resolved which is relevant for smectic A phases [601 ■ 

Second, the generalization to dynamics is in principle 
straightforward following the lines given in two spatial 



dimensions in Ref. [30[. For Brownian dynamics, the 
dynamical density functional theory |6ll - l63l | was gener- 
alized to orientational dynamics |4l| and can be used as 
a starting point to derive the order parameter dynamics. 
In general, the dynamics for -01 is conserved while that 
for uq and ■02 is not. However, though the derivation 
can be done in principle, it turns out in practice that 
the actual equations of motions for the order parame- 
ters involve a huge number of terms such that it is too 
tedious to state them all explicitly. The dynamical ex- 
tension could in principle be applied to the dynamics of 
topological defects Q and to interfacial dynamics near 
three-phase coexistence [65|- Also, the dynamics in ori- 
entational glasses [6a] could be explored. 

Third, polar particles with a head-tail asymmetry will 
violate the symmetry conditions used here for the direc- 
tor field. There should not be a principle obstacle to 
derive the free energy functional for polar particles as 
well. Furthermore, throughout the paper, we assumed 
uniaxiality. Biaxiality in the orientational distribution 
could also be included in the gradient expansions at the 
expense of more coefficients entering the picture. It will 
also be most interesting to see how the treatment of bond- 
orientational order in the framework of the PFC-model 
compares with the results available from continuum-type 
approaches. 

Finally it would be interesting to generalize the anal- 
ysis to active particles which are driven by a constant 
propagation speed along their orientation [63, ll^- A 
dynamical density functional approach was recently |69j 
proposed for active particles which could be used as a mi- 
croscopic starting point to derive systematically gradient 
expansions. 
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